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For a distribution advected by a simple chaotic map with diffusion, the "strange 
", eigenmode" is investigated from the Lagrangian (material) viewpoint and compared 

to its Eulerian (spatial) counterpart. The eigenmode embodies the balance between 
diffusion and exponential stretching by a chaotic flow. It is not strictly an eigenmode 
in Lagrangian coordinates, because its spectrum is rescaled exponentially rapidly. 
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There are two main types of coordinates used to represent fluid flow and 
dynamical systems. Eulerian (or spatial) coordinates are fixed in space, while 
Lagrangian (or material) coordinates follow parcels of fluid. Strange eigenmodes 
are persistent patterns in mixing — they can decay slowly, and hence remain vis- 
ible in the concentration field for a long time. So far, these have been studied 
U: from the Eulerian viewpoint. Here we describe the nature of the strange eigen- 
mode in Lagrangian coordinates for a simple map. It is not a true eigenmode 
because its wavelength is continuously rescaled in time. 

<N ; I. INTRODUCTION 

> 

The enhanced mixing of a passive scalar is one of the most direct consequences of chaos: 
a flow whose trajectories exhibit sensitivity to initial conditions will lead to rapid mixing. 
There are powerful theories based on the distribution of Lyapunov exponents [1-3] that 
link the mixing rate of the passive scalar with the chaotic properties of the flow. It has 
been recently suggested, following earlier work of Pierrehumbert [4-10], that the mixing 
properties of the flow can often be elucidated only by solving a full eigenvalue problem for 
the advection-diffusion operator, in an analogous manner to what is done for the kinematic 
dynamo [11]. The resulting eigenfunctions have been dubbed strange eigenmodes by Pier- 
rehumbert, and are closely related to Pollicott-Ruelle resonances in ergodic theory [12-14], 
which describe the long-time decay of correlations in mixing hyperbolic dynamical systems. 
Strange eigenmodes have also been observed experimentally [15, 16]. They are often called 
persistent patterns or large-scale eigenfunctions. 

The strange eigenmodes reflect a balance between advection and diffusion. On its own, 
advection is incapable of achieving mixing: it shuffles the concentration field of the passive 
scalar but does not decrease its fluctuations. The role of advection is to stir the concentration 
field, thereby creating sharp gradients in concentration. Physically, these gradients are 
reflected in the filamentation experienced by a blob of dye when it is stirred. The sharp 
gradients enhance the role of diffusion tremendously, and this allows mixing to proceed. As 
the diffusivity of the scalar is made smaller, the scale at which this mixing occurs decreases, 
so the concentration field appears very rough. In the limit of arbitrarily small diffusivity, the 
concentration field is not smooth: it consists of a superposition of strange eigenmodes. The 
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dominant one among these eigenfunctions is called the strange eigenmode, although there 
may be several of comparable importance. 

In the present work we tie the two types of theories together (i.e., Lyapunov exponent- 
based and strange eigenmode) for a specific system, which was already studied in [8] from the 
strange eigenmode viewpoint. The strange eigenmode represents a fundamentally Eulerian 
(spatial) view of mixing, whereas the Lyapunov exponent view is Lagrangian (material), 
following as it does the stretching history of fluid elements. In the strongly chaotic systems 
we deal with in the present context, Lagrangian and Eulerian coordinates are related by 
a convoluted transformation. This transformation is so complex that its specific form is 
inaccessible (even numerically) after some time, a reflection of sensitivity to initial conditions. 
For long times, the two frames must be regarded as essentially independent: we cannot 
simply solve a problem in Eulerian coordinates and transform to Lagrangian coordinates (or 
vice- versa). Thus we believe it is worthwhile to take a chaotic system whose solution has 
already been obtained in Eulerian coordinates and solve it in Lagrangian coordinates. As 
indicated above, this links two views of mixing together, and in particular illustrates the 
nature of strange eigenmodes in Lagrangian coordinates. This also indicates the source of 
the breakdown of local theories. We will show that there exists a kind of Lagrangian strange 
eigenmode, which is not quite an eigenmode but which exhibits similar features: specifically, 
it is an eigenmode if an appropriate time-dependent rescaling of coordinates is performed 
(exponential in time). This rescaling is closely related to the "cone" involved in diffusive 
problems in the presence of an exponentially stretching flow [17]. We call it the cone of 
safety, because modes inside it are sheltered from diffusion at a given time. 

We introduce the system to be studied, the perturbed cat map, in Section II. We find 
its finite-time Lyapunov exponents and eigenvectors using first-order perturbation theory. 
In Section III we partially solve the advection-diffusion equation for our map, again using 
perturbation theory. Some numerical work is needed to complete the solution, and this is 
described in Section IV. Finally, a few concluding remarks are offered in Section V. 



The strange eigenmode has for the most part been studied in maps, because these present 
great advantages for analytical work. This is also reasonable since experimental work has 
so far focused on time-periodic flows [15, 16]. As in [8], we consider the map 



defined on the unit two-torus, T 2 = [0, l] 2 . Here M is a matrix of integers with unit deter- 
minant, and 4>(x) is a doubly-periodic function, so that M is a diffeomorphism; specifically, 
we take 



where <f>(x) is chosen such that M is area-preserving. For K = 0, (2.1) is the usual cat 
map of Arnold [18]. The map (2.1) inherits much of the simplicity of the cat map, but the 
perturbation allows for more complex — and less singular — behavior. The action of the map 
is depicted in Fig. 1. For small K, the map is very close to the cat map, but the implications 
of the perturbation for mixing are profound, as we will now discuss. 

We are interested in the mixing properties of the map M. In Ref. [8] mixing in this 
map was investigated from an Eulerian perspective: advection alternated with diffusion and 
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FIG. 1: Action of the perturbed cat map for K = 0.4. (a) Initial pattern; (b) First iterate; (c) 
Second iterate; (d) Third iterate. The gray shading shows the action of the unperturbed cat map. 

the central object was the distribution of the concentration field in Eulerian coordinates. 
For K = 0, mixing occurs superexponentially in the map, due to the lack of dispersion 
in Fourier space. The concentration in a given Fourier mode is mapped entirely to one 
mode of higher wavenumber, and so on to ever higher wavenumbers. This sequence of 
wavenumbers have exponentially-growing magnitudes. Because diffusion is exponential in 
the wavenumber, the net result is superexponential decay (i.e., the exponential of minus an 
exponential in time). 

For K ^ 0, the situation is radically different. The map now disperses concentration 
among many Fourier modes at each iteration. In particular, some concentration is always 
mapped back to the lowest allowable wavenumber (the grave mode), on which the weak 
diffusion is ineffective. The decay of the scalar is then limited by how much concentration is 
mapped to this grave mode at each iteration. The grave mode thus forms the seed of a strange 
eigenmode, since an eigenmode is by definition a recurring feature. The concentration field 
thus settles into the slowest-decaying eigenfunction of the advection-diffusion operator — the 
strange eigenmode — analogously to earlier work [4-7, 9, 10]. 

Here we wish to solve the same problem as in Ref. [8], but in Lagrangian coordinates. 
This means that we focus on following fluid elements and describing how they deform under 
the action of the map. To solve the advection-diffusion problem in Lagrangian coordinates, 
it is necessary to have expressions for the finite-time Lyapunov exponents of the map [19] 
(or equivalently the coefficients of expansion) and their associated characteristic directions, 
as a function of Lagrangian coordinates, and not just their distribution (we will see why this 
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is so in Section III). Because the finite-time Lyapunov exponents are easily derived for the 
cat map (K = 0), we shall proceed perturbatively, assuming K is small. 

First let us give the Lyapunov exponents and associated characteristic directions for the 
unperturbed cat map. The Lyapunov exponents are the logarithms of the eigenvalues of M, 
and the characteristic directions are the corresponding eigenvectors. It is convenient to 
introduce an angle 9 in terms of which the eigenvectors of M are 

/. fcosd -sin6>\ , ■ 

( U S ^={sin9 cos 9 ) ^ 

with cos = 1(1 + l/y/E) and sin 2 9 — |(1 — 1/y/E). Then the corresponding eigenvalues 
of M are 

A u = A = |(3 + y/b) = 1 + cot # = cot 2 9, (2.4a) 
A s = A" 1 = |(3- Vh) = 1 -tan0 = tan 2 #. (2.4b) 

These equalities are specific to this particular angle, as is the relation tan 9 = cot 9 — 1. 
The u direction is associated with stretching, and s with contraction. 

The coefficients of expansion (given by A n and A~ n after n iterations of the map) and 
characteristic directions for the linear cat map are uniform in space. Now we derive their 
value for K nonzero but small, using perturbation theory. The problem is to find the 
eigenvalues and eigenvectors of the matrix g( n \ with components 



An) ._ ST dx ™ dx ™ ( 2 <\ 

a Pq - 2^ dXn dx n ' 
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often called the metric tensor (or Cauchy-Green strain tensor in fluid mechanics). Here X is 
the Lagrangian label (coordinate) and x = x^ n \X) is the nth iterate of the point X under 
the action of the map (2.1) (so that x(°\X) = X). The metric tensor describes the stretch- 
ing experienced at the nth iteration by a fluid element initially at X. Its eigenvalues and 
eigenvectors give the shape and orientation at the nth iteration of an ellipsoid representing 
an initially spherical infinitesimal element of fluid. 

Before we can apply perturbation theory to the metric tensor, we must find the form of 
the perturbation itself. To first order in K, the nth iterate of the map (2.1) is 



n— 1 



x (n) = M «(x) = M n • X + M ™ ' </>(M n_m_1 X). (2.6) 



m=0 

The Jacobian matrix of this transformation is 

n-l 



— = M n + V M m • -^-(M^ m - 1 X) • M"-" 1 " 1 (2.7) 

dX ^ dx y ' v ' 



m=0 



where the numerator corresponds to rows and the denominator to columns of a matrix. We 
must now construct the metric tensor (2.5), which to leading order in K is 

gP = M n M n + {h {n) + h&} (2.8) 
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with 

h (n) := m™ • V M m • — ^(M n_m_1 X) • M"-™- 1 (2.9) 

ox 

m=0 

and the tilde denotes the transpose of a matrix. The unperturbed metric is = M n M n , 
with eigenvalues A 2 ™ and A~ 2n and eigenvectors given by (2.3). The perturbation is the 
bracketed term in (2.8). Finding the eigenvalues and eigenvectors of the symmetric ma- 
trix to leading order in K is a straightforward application of perturbation theory for 
symmetric matrices, familiar from quantum mechanics. For more details we refer the reader 
to standard texts on the subject [20, 21]. 

To leading order in K, the coefficient of stretching is written as 

AP(X) = A n (l + Krj^(X)) (2.10) 

where A n is the coefficient of stretching of the unperturbed cat map, and the correc- 
tion A n Kr]^ n \X) is obtained from the perturbation by contraction with the unperturbed 
eigenvectors [20, 21], 

A n K^ n) {X) = \ u • {h {n) + fcW} ■ u = u • h {n) • u. (2.11) 

Observe that because the coefficient of stretching is the square root of the largest eigenvalue 
of gjp, there is an extra factor of 1/2 to leading order in K. Using the fact that, for M 
symmetric, M-u = u- M = Au, we find Eqs. (2.9) and (2.11) give 

n-l 

= S in 6 cos 9 cos (27r(M m • X) ± ) (2.12) 

m=0 

where we have substituted the specific form of the map, given by (2.2). The subscript '1' 
in (2.12) indicates the X\ component of a vector. 
The perturbed eigenvectors can be written as 

uP(X) = u + KC (n \X) s , sP(X) = s - KC (n \X) u , (2.13) 

where may be regarded as a small angle of rotation. Again we follow standard matrix 

perturbation theory [20, 21], so that the angle of rotation is given by 

K&\X)= \ 2n _ A J , (2.14) 
which after some reduction and the use of (2.2) yields 

c ( " ) = xct^ (c + )+c - >) - (2 - 15) 

with 

n-l 

Ci n) := i(cos2^ T l)^A ±2(n - m) cos(27r(M m -X) 1 ). (2.16) 

m=0 
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FIG. 2: Numerical solution (solid) and first-order solution r/ n ) (dashed) from Eq. (2.12) 
for K = 10~ 5 . The solutions diverge after several iterations because we are perturbing off a chaotic 
trajectory. 

Note that the asymptotic direction (n >> 1) is dominated by (+ \ so that 

n-l 

C (n) ~cos 2 ^^A- 2m cos(27r(M m -X) 1 ), n > 1. (2.17) 

m=0 

In this form it is easy to check that u • = s • V?/ n ) for n 1, as required by the 

differential constraint V • + • V log = [22-24] . (The derivatives are taken with 
respect to the Lagrangian coordinates X.) 

The first-order perturbative solution for the coefficient of stretching r/ n ) is compared to 
numerical results in Fig. 2, and similarly for the eigenvector in Fig. 3. Unlike the 
perturbed coefficients of stretching, which eventually diverge from the numerical solution 
because of the sensitivity to initial conditions, the perturbed eigenvectors converge very 
rapidly and are always close to the numerical result. 

It will be more convenient to express the metric tensor in terms of its eigenvalues and 
eigenvectors. The metric tensor (2.5) can be written in terms of the coefficients of expansion 
and characteristic directions as 

g {n) = [ A M] V n) U W + [A(")]-2g(n)g(n) 

To leading order in K, we have 
gP = A 2n uu+A- 2n ss+2K^ n \A 2n uu-A- 2n ss)+K(^ (A 2n - A~ 2n ) (us+su). (2.19) 




(2.18) 
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FIG. 3: Numerical solution (solid) and first-order solution (dashed) for the first component of the 
asymptotic eigenvector with K = 0.1, using the asymptotic result (2.17) in (2.13). The error 
is of order K 2 . 

The only dependence on X in (2.19) is contained in rf n ^ and C}~ n \ 

III. ADVECTION AND DIFFUSION 

Having derived the coefficients of expansion and characteristic directions of stretching 
(to leading order in K), we can now solve the advection-diffusion equation in Lagrangian 
coordinates. We will first discuss the case for an incompressible flow, and then make the 
transition to a volume-preserving map. 

The advection-diffusion equation for the concentration of a scalar, Q(x,t), advected by 
an incompressible velocity field v is 



d t e + v-d x o = Kd 2 x e, 



(3.1) 



where k is the diffusion coefficient. We define the transformation x(X,t) from Lagrangian 
coordinates X to Eulerian coordinates x by 



x = v(x,t), x(X,0) = X, 



(3.2) 



where the overdot denotes a time derivative at fixed X. We can then transform (3.1) to 
Lagrangian coordinates X [19, 22], 



Q = d x - (B-d x Q), 



(3.3) 
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where we reused the same symbol for Q(X,t). The anisotropic, nonhomogeneous, time- 
dependent diffusion tensor D is 

0:=k, , (3.4) 

where (7 is the metric tensor encountered in Section II. By construction, the advection 
term has disappeared from (3.3). The flow enters (3.3) indirectly through the metric tensor 
in (3.4), reflecting the enhancement to diffusion due to the deformation of fluid elements [19, 
22]. 

We now make the leap from a flow to a map: because the velocity field does not enter (3.3) 
directly, we may regard the time dependence in D as given by a map rather than a flow, 
and use the metric (2.5) in the diffusion tensor D [27]. We also write 9 (n) (X) for Q(X,t), 
where n denotes the nth iterate of the map. 

Since our map is defined on the torus, we can expand Q^ n \X) in Fourier components 
the resulting map, obtained by first Fourier transforming and then solving (3.3), is 

©L n) = E ex p(s (n) )^ n_1) > (3-5) 

where 



g<£) = -4tt 2 T / (k • D (n) • I) e~ 2 ^ k - £ > x d 2 
Jr 2 



x, (3.6) 



with T the period of the map. This is an exact result, but the great difficulty lies in 
calculating the exponential of Q^ n \ Again, we shall accomplish this perturbatively. 
For the torus map introduced in Section II, from Eq. (2.19) we obtain 

[gP]- 1 = A 2n ss + A- 2n uu + 2if?/ n) (A 2n ss- A^uu) - K ( {n) (A 2n - A~ 2n ) (us + su), 

(3.7) 

to leading order in K, where the only functions of X are rf n ^ and Inserting (3.7) 

into (3.6), we find 

Sffi = 42 + KB™ (3.8) 

where 

4? := -e (A 2 ^ 2 + A- 2n k 2 u ) 5 ki (3.9) 

and 

B<$ := -e (2 (A 2n k s i s - A~ 2n k u Q rj™ - (k u £ s + k s i u ) ((^^ + C W fc ,)) ■ (3. 10) 
Here we have defined 

e:=4n 2 nT (3.11) 
to agree with the notation in Ref. [8], as well as 

k u := k ■ u, k s :— k • s, (3.12) 

and similarly for £ u and £ s . 

Upon making use of the Fourier-transformed (2.12), (2.15), and (2.16) in (3.10), we find 

n-l 

Bkl = ~\ e ^ ^fe^ (<W+ei-M™ + ^fc,£-ei-M m ) (3.13) 
m=0 
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where ei is a unit vector in the Xi direction, and 

ft™ := sin 26 (A 2n A; s £ s - A~ 2n k n £ u ) + (k u £ s + k s Q (A 2 («" m ) s i n 2 6 - A~ 2 ( n ~ m ) cos 2 9) . 

(3.14) 

To obtain the full solution, we must now exponentiate (3.8) to give the transfer matrix 
in (3.5). Fortunately, for A diagonal there is a simple expansion, 

[exp(AH + KBV>)\ M = S M + KE k f, := B$$ ^ ~ ^ , (3.15) 

A kk ~ A ut 

valid to first order in K. We say fortunately because without such a formula it is very 
difficult to compute this matrix exponential — even numerically — due to the large size of the 
matrices (i.e., infinite) and their magnitude (i.e., growing exponentially in time). 

The A 2n term in AjjJ seems to imply that 0^ decays superexponentially fast as 
exp(— e A 2n kl). From Eulerian considerations [8], we know that for K ^ the decay is 
actually exponential after a short superexponential transient. This is because the term 
must be taken into account: it breaks the diagonality of 9^ n \ so that given some initial set 
of wavevectors, the concentration contained in those modes can be transferred elsewhere. In 
particular, it can transfer concentration to modes aligned with the unstable direction. We 
will see how this avoids superexponential decay in Section IV. 



IV. NUMERICAL RESULTS 
A. The Numerical Method 

At this point, solving (3.5) and (3.15) numerically seems like the only way forward. 
Clearly, attempting the solve this on a grid in Fourier space is hopeless: very high wavenum- 
ber modes are quickly populated so the resolution is exhausted very rapidly. Instead, the 
procedure we use involves keeping track of a list of excited Fourier modes (i.e., those that 
are nonzero to machine precision). We now describe this scheme. 

First, an initial wavenumber is seeded with some initial concentration. This mode will 
be damped by the diagonal part of the matrix in (3.15), and will also excite two new 
modes as seen in (3.13). Repeating this, starting now from three modes, we see that the 
number of excited modes grows exponentially. Thus it would seem that this procedure 
is not very advantageous; however, after a few iteration the diffusivity (the diagonal part 
in (3.15)) will damp most modes because is growing exponentially. Thus the modes 
that have been damped beyond redemption can be removed from the list. In this manner the 
number of excited modes eventually reaches a constant, though they consist of ever higher 
wavenumbers. Thus one can think of a "packet" of modes cascading through Fourier space 
towards larger wavenumbers. It is this packet that is the Lagrangian analogue of the strange 
eigenmode in Eulerian space, as we will discuss in Section IV B. 

Let us first present some results. Figure 4 shows the decay of the scalar variance 
for K = 0.001 and four values of the diffusivity. The agreement with the Eulerian re- 
sults is excellent for early times, but inevitably breaks down later. (In fact the variance 
eventually begins to increase, which is forbidden.) The agreement is also worse for smaller 
diffusivity. Both of these disagreements are a manifestation of the wavenumber dependence 
of the perturbation in (3.15): for k too large the perturbation becomes large, invalidating 
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the approach. Nevertheless, Fig. 4 clearly validates the calculation for times that are not 
too long. 

Another validation is shown in Fig. 5, where we plot the difference between the Eulerian 
and Lagrangian results as a function of K. The difference clearly scales as K~ 2 , showing 
that the two agree at leading order, as required for a first-order asymptotic result. 

We now interpret our results in greater detail, and look for a manifestation of the strange 
eigenmode in Lagrangian coordinates. 



B. The Lagrangian Strange Eigenmode 

The mechanism described in Section IV A is similar to that originally introduced (in the 
context of the kinematic dynamo problem) by Zeldovich et al. [17]: they basically solved the 
advection-diffusion equation in Lagrangian coordinates for a linear velocity field, and found 
that in order to avoid rapid superexponential decay one had to restrict attention to a "cone" 
of wavenumbers that are closely aligned with the unstable manifold of the flow (a similar 
approach was used later in Refs. [1-3]). The exponential shrinking in time of this "cone of 
safety" leads to an exponential decay of scalar variance at a rate given by the Lyapunov 
exponents. 

The problem with that approach is that a linear velocity field offers no possibility of 
dispersion in Fourier space. The wavenumbers in the cone of safety must have some con- 
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FIG. 5: Relative difference in the variance between the Eulerian and perturbative Lagrangian result 
as a function of K, at the fourth iteration. The dashed line indicates an K 2 dependence, showing 
that the two agree to first order in K. 

centration associated with them initially. What our numerical results show is that if one 
considers dispersion in Fourier space (of the type allowed by the term in (3.15)) then 
it is possible for concentration to be moved inside the cone from elsewhere. The Lagrangian 
equivalent of the strange eigenmode must live within the cone of safety, otherwise it would 
decay away superexponentially. But unlike Ref. [17] the decay rate in the present case is 
not determined by the shrinking of the cone: it is set by how much variance gets transferred 
into the cone at each iteration. 

Figure 6 shows a plot of the power spectrum of concentration. The magnitude of the 
concentration (log 10 ) is plotted vs the magnitude of the wavenumber normalized by ||fe|| max 
(its maximum value), which is proportional to A n . The concentration is normalized at each 
iteration such that the mode with largest concentration has unit magnitude. The iterations 
plotted are n — 6, 7, 8, 9, 10, 11 (circles) and n = 12 (black dots). Most of the circles appear 
as large black dots, because all these points lie on top of each other. Hence, the concentration 
is in an eigenstate, given that the wavenumber has been rescaled by a factor proportional 
to A n (i.e., such that the dominant peak is at unit rescaled wavenumber). This is what 
we interpret as the Lagrangian equivalent of the strange eigenmode (a good name might be 
"stretched eigenmode" in Fourier space). 

Some points in Fig. 6 exhibit a decay with iteration number (appearing as columns of 
circles, with higher iteration numbers lower on the graph): they belong to a more rapidly- 
decaying eigenfunction. Note that the peaks do not sharpen with iteration number, but 
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FIG. 6: The power spectrum of the Lagrangian strange eigenmode for n = 6, ...,11 (circles) 
and n = 12 (black dots). The large black dots are points that are the same for all iterations (after 
rescaling): this is the dominant strange eigenmode in Lagrangian coordinates. Both axes have 
been rescaled such that the dominant peak has unit amplitude and wavenumber (feu™ ^ is defined 
in Eq. (Al)). The hollow circles are due to an admixture of another, faster-decaying eigenfunction. 



more points are added to some of the tails. The eigenfunction appears extremely rough and 
discontinuous, though the peaks are indicative of some underlying continuum behavior. The 
seemingly isolated points actually tend to line up with a peak far below. Finally, note that 
the relative height (but not position or shape) of the peaks depends on K: the whole shape 
is stretched vertically as K is made smaller. This is because the term proportional to K 
controls the transfer of concentration "vertically" (with respect to Fig. 6) in the eigenmode 
at each iteration. 



V. DISCUSSION 

The Lagrangian strange eigenmode has some intriguing features: (i) It is rescaled expo- 
nentially in time, in order to remain within the cone of safety (so it is not a true eigenmode); 

(ii) Its power spectrum is very discontinuous, in sharp contrast to its Eulerian counterpart [8]; 

(iii) Its decay rate is set by how much concentration is moved into the "new" cone of safety 
at each iteration (since the cone is shrinking). In the Appendix we present an analytic result 
for a two-mode system which gives a simplified representation of the cone of safety. 

In the map analyzed here the exponential time-rescaling gives a proper eigenmode, since 
only the constant stretching is important at leading order. In a generic map the stretching 
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is a strong function of space, so the necessary time-rescaling would be position-dependent. 

It is hard to see how the Lagrangian approach presented here could be used in more 
realistic problems: perturbation theory was used extensively (which would not be applicable 
in most real situations), both for computing the finite-time Lyapunov exponents and the 
matrix exponential (3.15). We believe the approach is instructive nonetheless, giving as it 
does a picture of the strange eigenmode in Lagrangian coordinates. 

Our approach does not yield much information about the long-time behavior of the decay. 
There is currently a debate as to whether the mechanism presented in Refs. [1-3] gives a 
lower bound on the decay rate [25, 26]. Our perturbation expansion breaks down before this 
question can be answered. 
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APPENDIX A: THE TWO-MODE SOLUTION 

Though we have not found a general method of solution of (3.5) with the exponential 
given by (3.15), there is at least an approximate solution available that illustrates the broad 
features of a full solution. It also shows how the decay rate of the variance can become 
independent of the diffusivity in the Lagrangian viewpoint, as in Ref. [17]. 

The method is based on defining a class of "aligned" wavenumbers (i.e., that live inside 
the cone of safety), and retaining only two of these modes. These wavenumbers k^ are 
defined by 

k (p) = A p sin , k {p) = A- p cos 9, (Al) 

that is, fc (p) = M p • feo, where k o is any initial wavenumber for large enough p. Then k 
satisfies 

fc(p) _ fc(p-i) = e x ■ M p -\ (A2) 

so that with the choice k = k^ p \ £ = k^ 1 ^, the first Kronecker delta in (3.13) is unity 
for m = p — 1 . 

At the nth iteration, assume that only two wavenumbers are important: fc( n_d ) 
and fe( n_a!_1 ). The number d is a "lag" from the current iteration and will be adjusted 
later. 



15 



Define 

A d := exp 4t_ d)fe( „_ d) = exp (-e (A 2d cos 2 9 + A~ 2d sin 2 9)) (A3) 

which is independent of n, since we have defined p relative to the current iteration of the 
map. It can be shown that then the coupling from fet™-^ 1 ) to k^^ takes the simple form 

£ d := Bj!.,^,, = -\K (A d - . (A4) 

To recapitulate: at the nth iteration, the mode is mapped to itself with coupling 

amplitude A d , and ki n - d - 1 ) [ s mapped to fc( n_d ) with amplitude Z d . It is easy to show that 
in this two-mode situation the decay rate is determined by the magnitude of 8, d . All that 
remains is to find d. 

The "lag" d is obtained by maximizing E d over d; d has to be large enough that A 2d 
overcomes the tiny difiusivity in (A3) — so the two A terms don't cancel in (A4) — but not 
so large that A d is damped. We are thus justified in approximating A d ~ exp(— e A 2d cos 2 9) 
(the other term in (A3) is smaller by a factor A~ 4d , which is small even for d = 1). We then 
have 

£, = i^(T-T A2 ), T:=A*-i, (A5) 

since A d = (A d _i) A2 . This is easily extremized over T: the maximum \8. d \ is achieved 
for T = A~ 2 /( A _1 ) ~ 0.7198, for which \E d \ ~ 0.3074 K. The lag is then given by solving 
for d in terms of the extremizing T, 

d=l + - log ( ° g — ) ~ 0.5902 + 0.5195 log e~\ (A6) 

2 log A \ecos 2 9 J 

which scales logarithmically with the difiusivity Note that the decay rate is now completely 
independent of the actual value of the difiusivity: the lag adjusts itself to compensate, 
introducing a separation of scale between the dominant wavenumber fc*™ - ^ and the largest 
wavenumber in the system, k^ n \ 

The actual decay rate as e — > (from the Eulerian solution in [8]) is 0.5 K for small K, 
compared to the two-mode Lagrangian solution 0.3074 K. Thus, most of the important 
behavior is captured by the two-mode solution. The two modes can be seen in the spectrum 
of the strange eigenmode in Fig. 6: the dominant peak at ||fc||//ci n ^ = 1 is fc*™"^ (d — 2 
in this case), and the peak at ||fc||/A;i n 3 -* = A -1 ~ 0.3820 is k ( - n ~ d ~ 1 \ The other peaks are 
modes that could be included to get a more accurate expression for the decay rate. 

The two-mode solution also nicely illustrates the idea of the cone of safety: both modes 
are always inside it, and because the cone is shrinking by a factor A -1 at each iteration then 
the modes have to follow suit. The key difference with [17] is that here the concentration 
in the modes is mapped from one cone to another at each iteration, and is not part of the 
initial condition. 



